Search for: All records

Abstract: Atacama Large Millimeter/submillimeter Array observations have shown that candidate “post-starburst” galaxies (PSBs) at z~0.6 can retain significant molecular gas reservoirs. These results would imply that—unlike many model predictions—galaxies can shut down their star formation before their cold gas reservoirs are depleted. However, these studies inferred star formation rates (SFRs) either from [OII] line fluxes or from spectral energy distribution (SED) modeling and could have missed large dust-obscured contributions to the SFRs. In this study, we present Keck/NIRES observations of 13 massive (M_* >= 10^11M_⊙) PSBs, which allow us to estimate Hα SFRs in these gas-rich PSBs. We confirm the previously inferred low SFRs for the majority of the sample: 11/13 targets show clear Hα absorption, with minimal infilling indicating dust-corrected SFRs of <4.1Msun/yr. These SFRs are notably low given the large H2 reservoirs (∼(1–5) × 10^10Msun) present in 5/13 of these galaxies, placing them significantly offset from star-forming galaxies on the Kennicutt–Schmidt relation for star-forming galaxies. The [NII]/Hα ratios of all 13 PSBs imply contributions from non-star-forming ionization mechanisms (e.g., active galactic nuclei, shocks, or hot evolved stars) to their Hα emission, suggesting that even these low ongoing SFRs may be overestimated. These low Hα SFRs, dust corrected using Av estimates from SED fitting, confirm that these galaxies are very likely quiescent and, thus, that galaxies can quench before their cold gas reservoirs are fully depleted. 

Abstract We explore the cosmic evolution of the fraction of dust-obscured star formation predicted by thesimbacosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up toz= 3, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, atz> 3 we find strong evolution at fixed galaxy stellar mass toward greater amounts of obscured star formation, in tension with high-redshift observations. We explain the trend of increasing obscuration at higher redshifts by evolving star-dust geometry, as the dust-to-stellar mass ratios remain relatively constant across cosmic time. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust-to-stellar mass ratios at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust-obscured star formation history peaks aroundz∼ 2−3, and becomes subdominant atz≳ 5. The dominance of obscured star formation at redshiftsz≲ 4 is consistent with our results for the evolution of the obscured star formation fraction at fixed stellar mass to higher values at higher redshift because there exist fewer massive, heavily obscured galaxies at high redshift. 

Abstract We presentslick(the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [Ci], and [Cii] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer codedespotic,slickcomputes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their physical properties and individual environments. We validate our results by applyingslickto the high-resolution run of theSimbasimulations, testing the derived luminosities against empirical and theoretical/analytic relations. To simulate the line emission from a universe of emitting clouds, we have incorporated random forest machine learning (ML) methods into our approach, allowing us to predict cosmologically evolving properties of CO, [Ci], and [Cii] emission from galaxies such as luminosity functions. We tested this model in 100,000 gas particles, and 2500 galaxies, reaching an average accuracy of ∼99.8% for all lines. Finally, we present the first model light cones created with realistic and ML-predicted CO, [Ci], and [Cii] luminosities in cosmological hydrodynamical simulations, fromz= 0 toz= 10. 

Abstract In the current era of high-zgalaxy discovery with JWST and the Atacama Large Millimeter/submillimeter Array, our ability to study the stellar populations and interstellar medium conditions in a diverse range of galaxies at Cosmic Dawn has rapidly improved. At the same time, the need to understand the current limitations in modeling galaxy formation processes and physical properties in order to interpret these observations is critical. Here, we study the challenges in modeling galaxy dust temperatures, both in the context of forward modeling galaxy spectral properties from a hydrodynamical simulation and via backwards modeling galaxy physical properties from mock observations of far-infrared dust emission. Using thesimbamodel for galaxy formation combined withpowderdayradiative transfer, we can accurately predict the evolution of dust at high redshift, though several aspects of the model are essentially free parameters (dust composition, subresolution dust in star-forming regions) that dull the predictive power of the model dust temperature distributions. We also highlight the uncertainties in the backwards modeling methods, where we find the commonly used models and assumptions to fit far-infrared spectral energy distributions and infer dust temperatures (e.g., single temperature, optically thin modified blackbody) largely fail to capture the complexity of high-zdusty galaxies. We caution that conclusions inferred from both simulations—limited by resolution and post-processing techniques—and observations—limited by sparse data and simplistic model parameterizations—are susceptible to unique and nuanced uncertainties that can limit the usefulness of current high-zdust measurements. 

Abstract This paper gives an overview of Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star formation (TEMPLATES), a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies, two extremely dusty and two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 ≤z≤ 4.2 galaxies, in particular Hα, Paschenα, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [Oiii] 5007 Å and Hβ; the other two targets have JWST coverage of polycyclic aromatic hydrocarbon 3.3μm and complementary Atacama Large Millimeter/submillimeter Array data covering the [Cii] 158μm emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has the technical goal to establish best practices for the integral field units within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES’s observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks. 

We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain-size distributions, and ionization states with a new model for dust evolution in galaxy simulations. We apply these models to three idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and a starburst disk. Our main results are as follows. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to shatter grains efficiently, driving up the fraction of ultrasmall grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain-size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to investigate fully the complex interplay of processes that drive PAH band luminosities in galaxies. 

Abstract One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the simba cosmological simulation with the prospector SED fitting code. We consider the impact of the commonly assumed uniform screen model and introduce a new nonuniform screen model parameterized by the fraction of unobscured stellar light. This nonuniform screen model allows for a nonzero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a nonuniform screen model, with the median offset in UV attenuation decreasing from −0.30 dex with a uniform screen model to −0.17 dex with the nonuniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the nonuniform screen model, decreasing the SFR offset on average by 0.12 dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.